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The phase diagram of a quantum XY spin chain with Gaussian-distributed random anisotropies 
and transverse fields is investigated, with focus on the fidelity susceptibility, a recently introduced 
quantum information theoretical measure. Monitoring the finite-size scaling of the probability dis- 
tribution of this quantity as well as its average and typical values, we detect a disorder-induced 
disappearance of criticality and the emergence of Griffiths phases in this model. It is found that 
the fidelity susceptibility is not self- averaging near the disorder-free quantum critical lines. At the 
Ising critical point the fidelity susceptibility scales as a disorder-strength independent stretched 
exponential of the system size, in contrast with the quadratic scaling at the corresponding point 
in the disorder- free XY chain. Along the line where the average anisotropy vanishes the fidelity 
susceptibility appears to scale extensively, whereas in the disorder-free case this point is quantum 
critical with quadratic finite-size scaling. 



I. INTRODUCTION 

In the last few years, tools from the field of quantum in- 
formation theory have found extensive use in the study of 
the phase diagrams of quantum systems. One such tech- 
nique, the fidelity approach to Quantum Phase Transi- 
tions (QPTs) has been successfully applied to various sys- 
tems possessing quantum critical points (see @ 
for a review) . This technique can be generalized to finite- 
temperature systems [gI, ^vj , classical phase transitio ns [p| , 
and topological phase transitions [§, [13, [ll|, [H, [H, [ij 

In a recent letter [l^ we have studied the fidelity in 
the context of disordered quantum systems. The physics 
peculiar to disordered quantum systems is reflected in 
the properties of the fidelity, a quantity not previously 
used to investigate such systems. Here we study the 
scaling behavior and provide details concerning the zero- 
temperature phase diagram of the disordered quantum 
XY model in a transverse field, a prototypical model in 
the context of disordered quantum systems. 

The paper is organized as follows. Section[n]is devoted 
to defining the model, along with a review of known re- 
sults about its phase diagram and the basics of the fi- 
delity approach. Section IIIII presents the numerical re- 
sults of our study and discusses the main features of the 
fidelity for disordered quantum chains. Our conclusions 
are presented in Section IIVI 



II. METHOD AND MODEL 

It is known that disorder can have interesting effects 
on a system's phase diagram pj6| . In particular, Griffiths 
phases may arise as a result of the randomness 17]. Here 
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we study the disordered anisotropic quantum XY spin 
chain in a random transverse field, a model where the 
disorder-free case can be analytically solved [3 and for 
which some exact results are known in the disordered 
case [l9|, |20|] ■ Its Hamiltonian is given by 

= - 2^ 2 ' 2 ' ' ' ^ ' 

i=l 

where aj^'^'^^ are Pauli matrices, and the fields 
and anisotropies 7; are independent Gaussian-distributed 
random variables. The average field and anisotropy are 
denoted by A and 7, respectively. The variance is taken 
to be the same for both the field and anisotropy distri- 
butions. 

The Jordan- Wigner transformation maps this system 
onto quasi- free spinless fermions [l^. Neglecting the 
boundary term and taking the system to be closed in 
the fermion index, we obtain a Hamiltonian of the form 

L ^ L 

^=J2 ^I^^JCJ + 2 ^ (4-B,j4 + C.Sy C,;) , (2) 

where A and B are symmetric and antisymmetric real 
Lx L matrices, respectively. Explicitly: Aij = —{2XiSij + 
Si,j+i + Si+i^j), Ail = Ali = -1 and = 'IjSij+i - 
7i(5i+ij, BiL = 7l = -Bli- 

The Hamiltonian may be rewritten in terms of the ma- 
trix Z = A — B, which contains all information about the 
system. Performing the polar decomposition of Z we ob- 
tain the matrices A and T such that Z — AT, where A 
is a positive semi-definite matrix and T is unitary. From 
the eigenvalues of A one obtains the single-particle energy 
spectrum [21*1. 

For systems at zero temperature, the fidelity is simply 
the absolute value of the overlap between ground states 
corresponding to nearby points in parameter space. Near 
a quantum critical point the ground state changes rapidly 
for small shifts in the tuning parameters, an effect which 
is reflected in a corresponding decrease of the fidelity. 
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The ground state fidelity can be cast in terms of the 
unitary matrix T in the following way [2^ 



F{Z,Z) = \\dct 



T + T, 



(3) 



where T and T are respectively the unitary parts of the 
matrices Z = Z{x) and Z = Z{x'), evaluated at the 
model parameters x and x' . The corresponding fidelity 
susceptibility is defined as [H, H^l 



lim ■ 



-2liiF{x,x + Ax) 



(4) 



and can be written in terms of the unitary matrix T as 

xi^) = l\\dxT\\l, (5) 

where || • \\f is the Frobenius norm. For a derivation of 
Eq. ^ see the Appendix. 

We evaluate the fidelity susceptibility using ([5]) by per- 
forming a singular value decomposition of Z 

Z = f/SF^ = {UY.U^){UV^) = AT, 

where U and V are unitary matrices. Note that the 
fidelity susceptibility is defined for infinitesimally sep- 
arated points along any chosen direction in parameter 
space. 



A. Disorder-free case 

Before considering the effects of disorder in the XY 
chain, let us first recall the behavior of the fidelity sus- 
ceptibility for the disorder-free case, where = A, 7^ = 
7,Vi e {1,...,L} [ij. The system can then be found in 
one of three phases. For |A| > 1 it is paramagnetic, and 
for |A| < 1 and 7 > (7 < 0) the system is ferromagnetic 
along the x- direction (j/-direction) . The boundary be- 
tween any two of these phases is a quantum-critical line 
corresponding to a second-order quantum phase transi- 
tion. Here, we refer to the transition driven by the mag- 
netic field as the Ising transition, and to the transition 
driven by the anisotropy coupling as the anisotropy tran- 
sition. At the quantum-critical points there is an avoided 
level crossing between the ground state and the first ex- 
cited state. As shown in Figs. [Tfa) and[Hb), one observes 
a maximum of the fidelity susceptibility at both the Ising 
and anisotropy critical lines. 

Moreover, the finite-size scaling dimension of the fi- 
delity susceptiblity in Figs. [Ijc) and[ljd) show it to be 
extensive away from criticality and superextensive (scal- 
ing quadratically with L) at the critical points. This 
scaling behavior holds for both the Ising and anisotropy 
critical lines. The apparent subextensive scaling in the 
immediate vicinity of the anisotropy critical point is a 
numerical artifact due to the narrowing of the fidelity 
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FIG. 1: (color online) The disorder-free case: (a) Fidelity 
susceptibility near the Ising transition, with 7 = 1 and system 
size L = 500. (b) Near the anisotropy transition, with A — 0.5 
and system size L — 500. (c) Finite-size scaling dimension of 
the fidelity susceptibility near the Ising transition, with 7 = 1. 
The sizes considered range from L = 100 to 700, with x ~ 
L^^ . (d) Finite-size scaling of the fidelity susceptibility near 
the anisotropy transition, with A = 0.5. The sizes considered 
range from L = 100 to 600. 



susceptibility peak as the system size grows. The Ising 
transition does not show this behavior, since the narrow- 
ing appears to occur more slowly than for the anisotropy 
transition. 

It has been shown [23] , for translationally invariant sys- 
tems, that superextensive finite-size scaling of the fidelity 
susceptibility implies a vanishing gap and therefore quan- 
tum criticality. As a result, for clean systems the points 
in parameter space corresponding to superextensive scal- 
ing of this quantity mark quantum critical regions. When 
randomness is introduced, translational invariance is lost, 
and hence superextensive scaling of the fidelity suscepti- 
bility does not necessarily imply quantum criticality. As 
discussed before, we find that locations of superexten- 
sive scaling reveal more general behavior beyond quan- 
tum criticality, namely Griffiths phenomena [17| . 



B. Random XY chain 

The effects of disorder on the physics of quantum mag- 
nets has been studied mainly using the strong-disorder 
renormalization group technique (SDRG) [25l. [26l[27j . A 
different approach has been used in the work of McKen- 
zie and Bunder [l^, , where the critical behavior of the 
disordered XY chain has been studied using a mapping 
to random-mass Dirac equations. The properties of the 
solutions of these equations imply the disappearance of 
the anisotropy transition in the presence of disorder. Fur- 
thermore, Griffiths phases are predicted to appear both 
around the Ising critical line and the anisotropy 7 = 
line. These results, together with the analysis performed 
by Fisher [2^, [23| , are significant since they analytically 
show the drastic effects that disorder can have on the 
critical properties of a quantum system. 
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At fixed 7 the XY random chain is closely related to 
the random trans verse- field Ising chain (RTFIC), which 
is anoth er p rototypical model for disordered quantum 
systems [27|- Since the RTFIC is representative of the 
universality class of Ising transitions for all values of 7, 
let us review what is known for this model. The Hamil- 
tonian of the RTFIC is H ^ - J2f^„^ [J^crfcrf+l + h^af] , 
where Ji and hi are random couplings and fields respec- 
tively. The system is critical when the average value of 
the field equals the average value of the coupling. Us- 
ing the SDRG one obtains that, at the quantum critical 
point, the time scale r and the length scale L are related 
by Inr ~ L^/^. This results in an infinite value for the 
dynamical exponent z at criticality [27l | . The distribution 
of the logarithm of the energy gap e at criticality broad- 
ens with increasing system size, in accordance with the 
scaling relation Ine ~ —L^/'^ [2g|. In the vicinity of the 
critical point the distribution of relaxation times is broad 
due to the presence of a Grifhths phase, characterized by 
a non-universal dynamical exponent z depending on the 
distance from the critical point. This dependence can be 
used as an indicator for the Grifhths phase. 

In [1^ a study of the phase diagram of a random XY 
spin chain in a random transverse field was performed 
using the fidelity approach. There it was shown that 
superextensive finite-size scaling of x signals the presence 
of a quantum phase transition close to the Ising critical 
line, while the minimum in x close to the anisotropy 7 — 
line is consistent with the absence of a phase transition 
in that parameter region, see Figs. Ufa) and (b). 



nent z and the scaling dimension of y [24 



A, 



2z + 2 - 2A, 



o 



(6) 
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FIG. 2: (color online) (a) Average fidelity susceptibility near 
the Ising transition, for L = 500, 7=1, and o — 0.3, (b) 
Average x near the 7 = line, for L = 500, A = 0.5, and a = 
0.3, (c) Scaling dimension of x for the same set of parameters 
near the Ising line, considering L G {100,200,300,400,500}, 
(d) Scaling dimension of x for the same set of parameters near 
the 7 = line, considering L e {100, 200, 300, 400, 500}. 



The Griffiths phases of the model manifest themselves 
in a non-universal dependence of the finite-size scaling di- 
mension of the fidelity susceptibility on the distance 
from the disorder- free critical point, see Figs. HJc) and 
(d). The relation between the dynamical scaling expo- 



establishes the connection between the fidelity suscepti- 
bility and the Griffiths phase. In Eq.® Aq is the scaling 
dimension of the relevant operator driving the transition. 
Eq. ([6]) implies a non-universal scaling dimension for 
the fidelity susceptibility A^, provided that the behavior 
of the unknown scaling dimension Aq does not exactly 
cancel that of the dynamical exponent. In the following 
we would like to study, using other methods, the extent 
of the Griffiths phase for this model. In particular, we 
would like to verify that the entrance into a Griffiths 
phase is indeed reflected by a changing scaling behavior 
of the fidelity susceptibility. 

In our numerical analysis we consider system sizes 
L £ (100, 200, 300, 400, 500} and for each system size we 
compute 10^ disorder realizations. We take the exter- 
nal fields (anisotropics) to be independent and identically 
distributed Gaussian random variables with standard de- 
viation a and mean A (7) . We consider the range of val- 
ues (0.1, 0.2, 0.3, 0.4} for the standard deviation cr. This 
disorder strength can be considered strong with respect 
to the value of the other parameters. We denote with 
Have ^^"^ arithmetic mean over all 10^ disorder realiza- 
tions. 

The width of the Griffiths phase for the XY model with 
weak Gaussian disorder in the continuum limit is known 
due to the work of McKenzie [l^. Let us denote the 
distance from criticality with S, where ^ = corresponds 
to the points in parameter space where the pure system is 
critical. From 131 one can compute that near the Ising 
transition, where the field A drives the transition, 6 = 

^-^^^^T-^, while near the anisotropy line 5 — ^ ■ 

For the so-called commensurate case, which corre- 
sponds to the Ising transition for this system since we 
have disorder in both the field and anisotropy, McKen- 
zie showed that the disorder-averaged density of states 
[p(-^)]avc //^o diverges at zero energy within the range 
1^1 < 1/2 away from the critical point [l^. Here po is the 
"high-energy" density of states. This divergence implies 
that the gap distribution function P{AE) also diverges 
at zero energy, since a large density of states means a 
vanishingly small gap. Note that a diverging probability 
of having a vanishing gap does not necessarily imply that 
the density of states is also divergent, since the gap dis- 
tribution only provides information about the position 
of the first excited energy level relative to the ground 
state energy. However, since a diverging gap distribution 
should be expected to occur as a result of a divergence 
of the low-energy density of states, we will use it to give 
a rough estimate of the extent of the Griffiths phase. 

For the incommensurate case, which includes the 
anisotropy transition, [p{E — O)]^^^^, /po is of the order of 
unity for some range of parameters about 5 — Q, imply- 
ing effective gaplessness, and [p{E = 0)]^^^ /po is much 
smaller than unity for \5\ >> 1, giving an effectively fi- 
nite gap [1^. Note that the boundary of the Griffiths 
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phase in this region is expected to be less defined than 
near the Ising transition, since the zero-energy density of 
states does not diverge at any value of 6. 

These results apply for the case of weak disorder, 
but we consider a range of moderate to strong disor- 
der strengths. In order to compare with the results of 
McKenzie for the Griffiths phase extent, we propose a 
rough criterion for determining the extent of the Grif- 
fiths phase using the gap distribution. Assume that the 
Griffiths phase lies within the range of parameter val- 
ues for which the distribution P{AE) has a maximum 
for AE = 0. At some point away from criticality the 
distribution maximum moves away from zero, eventually 
becoming approximately Gaussian far from the disorder- 
free critical point. We have determined the range of pa- 
rameters for which P{AE) has a maximum at zero gap. 
The width of this range of parameters is independent of 
system size and scales with the variance of the coupling 
or field distributions, in accordance with the result of 
p^ . However, this criterion gives an extent several times 
larger than the McKenzie value, around both the Ising 
and 7 = lines. This difference may be due to strong 
disorder or, rather, the result of an overestimation of the 
Griffiths phase extent since the distribution P{AE) hav- 
ing a maximum at zero does not directly imply a diver- 
gent zero-energy density of states. Nonetheless, consis- 
tency between different estimations of the Griffiths phase 
extent support the validity, at least qualitatively, of this 
approach. 

Fig-® shows the gap distribution in the vicinity of 
the 7 = line. For small anisotropics the distribution 
has a pronounced peak at zero gap. Moving away 
from the anisotropy line, the distribution develops 
a peak slightly away from AE = 0, and for larger 
anisotropics the distribution becomes Gaussian with a 
vanishing probability of having zero gap. Similar behav- 
ior holds for the Ising transition as the average field A 
is adjusted away from the finite-size pseudocritical point. 
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FIG. 3: (color online) Gap distribution near the anisotropy 
line, with A = 0.5, a = 0.3, and system size L — 500. Plotted 
distributions are for 7 = 0.14, 0.18, 0.22, and 0.30. 



III. RESULTS 
A. Average and typical values 

In computing the average of some physical quantity, 
namely the arithmetic mean over many disorder realiza- 
tions, any rare but large values will significantly affect 
the result. On the other hand, the geometric mean over 
disorder realizations gives a more representative measure 
of the typical values of the physical quantity. Recall that 
the arithmetic mean provides an upper bound for the geo- 
metric mean when averaging over a set of positive values, 
and the two are equal only when taking the mean of a 
constant set of values. 

To observe the presence of large fiuctuations in the 
fidelity susceptibility, we plot in Fig. [UJa) the disorder- 
averaged as well as typical fidelity susceptibility in the 
vicinity of the Ising critical point for a disorder strength 
a = 0.3. Notice that in the vicinity of the critical point 
the average becomes significantly larger than the typical 
value, indicating that there are instances of large fidelity 
susceptibilities that skew the arithmetic average towards 
a greater value. 

Near the anisotropy line, as shown in Fig. [Hb), there 
are also regions where the average fidelity susceptibility 
becomes much larger than the typical value, but now 
the positions of largest difference do not correspond to 
a critical point. Indeed, at the point 7 = 0, which in 
the disorder-free case is critical, the average-typical dif- 
ference is much smaller than it is at the two offset peaks. 
This is evidence for the disappearance of the anisotropy 
transition as a result of the disorder. 

For translationally invariant systems, which are 
disorder-free, it has been shown that superextensive 
finite-size scaling of the fidelity susceptibility implies 
quantum criticality The disordered XY chain 

does not have translational invariance, so locations of 
superextensive scaling do not necessarily imply criti- 
cality. However, it is still useful to consider finite-size 
scaling, since comparison with the disorder-free case 
may suggest in what way the phase diagram changes 
as a result of disorder. Fig.® shows the finite-size 
scaling dimension of the typical fidelity susceptibility 
near the Ising transition. The locations of the maxima 
of the typical fidelity susceptibility and finite-size scaling 
dimension coincide, and are shifted from the pure 
pseudocritical point due to finite-size effects. Also, the 
maximum scaling dimension obtained with disorder is 
smaller than the pure case of quadratic scaling in L, and 
this maximum value decreases with increased disorder 
strength. 

In Fig. ([5]) , the finite-size scaling dimension of x around 
the anisotropy line 7 = is shown for the same set of 
disorder strengths. Notice that the scaling depends on 
distance from the 7 = line and at 7 = the scaling 
is approximately extensive, as it is when far from the 
anisotropy line. For sufficiently small disorder and sys- 
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FIG. 4: (color online) Average (blue) and typical (red) fidelity 
susceptibility about: (a) the Ising line, for 7 = 1, L = 500, 
and a — 0.3 (b) the 7 = line, for A = 0.5, L = 500, and 
a = 0.3. 




FIG. 5: (color online) Finite-size scaling dimension of the 
typical fidelity susceptibility about the Ising transition, for 
7 = 1, for a £ {0.1, 0.2, 0.3, 0.4}. Scaling fit considered system 
sizes L = 200, 300, 400, 500, with x ~ assumed. 




FIG. 6: (color online) Finite-size scaling dimension of the 
typical fidelity susceptibility about the anisotropy line, for 
A = 0.5, for a £ {0.1,0.2,0.3,0.4}. Scaling fit considered 
system sizes L = 100, 200, 300, 400, 500, with x 
assumed. 




tem size there may appear to be only a single peak in the 
typical fideUty susceptibility at 7 = 0, suggesting that 
for these finite systems emergent criticality is still felt 
even though the quantum phase transition in the ther- 
modynamic limit disappears as a result of the disorder. 
However, increasing the system size reveals a double peak 
with a peak offset which grows with the strength of the 
disorder. In both the Ising and anisotropy regions, the 
width of the parameter interval giving scaling dimensions 
larger than a particular value scales approximately with 
the variance of the disorder distribution. This scaling 
behavior agrees with that given by the previously men- 
tioned gap distribution criterion for the Griffiths phase. 

In Fig. ([7]) the disorder-averaged gap is plotted for the 
four disorder strengths. The behavior of the gap corre- 
sponds closely with that of the fidelity susceptibility, in 
that the average gap minima have the same location as 
the typical x maxima. Note that the vicinity of the Ising 
critical line and the 7 = line are regions of effective 
gaplessness which we associate with quantum criticality 
or Griffiths phases. 

An experiment on such a disordered system would 
consider only one particular realization of disorder, and 



FIG. 7: (color online) (a) Disorder- averaged gap near the Ising 
transition for disorder strengths a £ {0, 0.1, 0.2, 0.3, 0.4}, L = 
500, and 7=1. (b) Disorder-averaged gap near the 7 = 
line for the same range of disorder strengths, L — 500, and 
A = 0.5. 



as a result would not necessarily observe the disorder- 
averaged value of an observable but rather a typical 
value. Here, we would like to study whether a measure- 
ment of the fidelity susceptibility for a large system coin- 
cides with the average value, and to do this we must see 
for what conditions x is a self-averaging quantity [2^ . 
Consider the quantity R^{x,L) = Var[x(x)]/[x(a;)]ave- 
We expect that R^ix, L) for fixed x will scale as a power 
law in the system size L, Ry_ ~ L^. If 6 — —1 then 
we say x is self-averaging, if 6 < then x is weakly 
self-averaging, and if 6 > then x is not self-averaging. 
In Fig. ([8]) we indicate the regions for which x is self- 
averaging, weakly self-averaging and non-self-averaging 
for various disorder strengths near the Ising transition as 
well as the 7 = line. 
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FIG. 8: (color online) Regions about the Ising transition (a) 
and anisotropy line (b) for which x is weakly self- averaging 
(blue) and non- self- averaging (orange), as a function of disor- 
der strength o. Outside these regions x is self- averaging. 
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B. Fidelity susceptibility distributions 

1. Near the Ising line 

Far from the Ising critical line, the distribution of the 
fidelity susceptibility is Gaussian, see Fig.©. However, 
in the vicinity of the Griffiths phase and the critical point 
the distribution is non-Gaussian, developing a slowly- 
decaying tail towards large fidelity susceptibilities, as 
shown in Fig. p^ . This tail reflects the presence of rare 
but large fidelity susceptibilities, and is expected to arise 
either in a Griffiths phase or in a quantum critical region. 
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FIG. 9: (color online) (a) Distribution of x f^^r from the Ising 
critical point, with A = 1.49 (blue). Overlaid is a Gaussian 
distribution function with the same mean and variance (red) . 
(b) Distribution of In x for the same parameters (blue) along 
with the corresponding distribution function of the logarithm 
of a Gaussian random variable. Here a = 0.1, L = 500 and 
7 = 1. This distribution is well described as Gaussian. 

Now we explore how the distribution of the logarithm 
of the fidelity susceptibility changes as the system size 
is varied, with all other parameters fixed. Far from the 
Ising transition the distribution of ln(x) narrows with 
increasing system size, as shown in Fig. [TlTa). The po- 
sition of the peak of the distribution remains fixed for 
a rescaling X ^ xl^ (see Fig. [TlT b)). but the width of 
the distribution decreases slightly. Choosing a more gen- 
eral scaling assumption ln(x) \n{x/ L") allows for 
an improved collapse of the distributions in this region, 
indicated in Fig. IllT c). The fit parameter a essentially 
translates the distribution, while the parameter (3 adjusts 



FIG. 10: (color online) (a) Distribution near the Ising critical 
point, with A = 1.03 (blue). Overlaid is a Gaussian distri- 
bution function with the same mean and variance (red), (b) 
Distribution of In x for the same parameters (blue) along with 
the corresponding distribution function of the logarithm of a 
Gaussian random variable. Here a — 0.1, L = 500, and 7 = 1. 
Clearly this distribution is not Gaussian, as can be seen by 
the much slower dropoff of the tail towards large fidelity sus- 
ceptibilities. 



the width. Such a scaling would imply a size dependence 
X ^ L" exp , where a, (3 > 0. However, for asymptot- 
ically large system sizes this would lead to subextensive 
scaling, and would thus not be expected to hold for all 
L. It appears that this apparent non-power-law scaling 
for the range of sizes we have considered may be due to 
finite-size effects, and we speculate that an assumption 
of extensive scaling would lead to an improved collapse 
for sufficiently large system sizes. 

At the Ising pseudocritical point the distribution 
P(lnx) broadens significantly with increasing system 
size, and a rescaling ln(x) L^^ ln(x) gives a good 
collapse for a value of the fit parameter (3 = 0.26 (see 
Fig. p^ ). This collapse and value of fit parameter holds 
for the pseudocritical points corresponding to all four 
disorder cases we have considered. A rescaling of this 
kind suggests that the fidelity susceptibility scales as a 
stretched exponential of the system size at the critical 
point rather than quadratically as in the pure XY chain. 
For the random transverse field Ising chain [2^, it is 
known that the energy gap vanishes as AE ~ exp — \/L at 
the Ising critical point. Recalling the alternative expres- 
sion for the fidelity susceptibihty x = X^n/o ^ ^'^^1^^ " a" ^ 
[1^, we expect that the first term in this series would 
dominate, and that the fidelity susceptibility might scale 
as X ~ 1/{AE)^. However, this crude argument appears 
not to be consistent with the scaling of the energy gap of 
the RTFIM, perhaps because of a lack of universality in 
the power of L in the stretched exponential. 



2. Near the anisotropy line 

Just like for the Ising case, far from the anisotropy line 
the distribution P{x) is well-approximated as a Gaus- 
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FIG. 11: (color online) (a) Distribution P(lnx)) at A = 1.39 
for L G {200, 300,400,500}. (b) Distribution P(ln(x/i)). (c) 
Rescaled distribution, ln(L~"x)- Here a — 0.2, 7 = 1, and 
the fit parameters are a = 0.62 and f3 = 0.5. Including this 
exponential correction improves the collapse of the curves. 
This scaling would suggest x ~ expL"''. 
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FIG. 12: (color online) (a) Distribution P(lnx)) at the pseu- 
docritical point A = 1.02 for L G {200,300,400,500}. (b) 
Distribution after a power-law rescaling ln(L~"x)) where 
a = 1.58. (c) Distribution of L~'^lnx for f3 = 0.26 and 
the same parameters. Here a = 0.2 and 7 = 1, though this 
collapse applies for all other disorder strengths considered. 



sian. Closer to the 7 = line the distribution looks 
much like the distribution of Inx(A) in the vicinity of the 
Ising line, becoming non-Gaussian with a slowly-decaying 
tail towards large values of x- Fig. (fT3|) shows the distri- 
bution P(\n{x)) for 7 = and noise strength a = 0.3, 
and Fig. fM)) shows the same quantity for 7 = 0.03, the 
value of average anisotropy coinciding with the peak in 
the typical value of x for that magnitude of disorder. 

Considering the point 7 = 0, as the system size 
increases the distribution P(lnx) does not change width, 
so a rescaling X ^ x/L gives a good collapse, see 
Fig. (fT5|) . This scaling also agrees with the extensive 
scaling of the average fidelity susceptibility at 7 = 0. 
Moving 7 away from this point in either direction, 
soon the distribution begins to shift superextensively, 
as shown in Fig. (fT6|) . For all values of 7 in this peak 
region, a rescaling of the form x ^ x/ gives a 
good collapse, where a is the fit value of the finite-size 



scaling dimension of the corresponding typical fidelity 
susceptibility. Continuing to move 7 away from the peak 
in the typical fidelity susceptibility, the distribution 
begins to narrow slightly as in the off-critical Ising case. 
However, a rescaling X ^ x/L appears to give a good 
collapse for 7 sufficiently large in magnitude. 
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FIG. 13: (color online) (a) Distribution of x at the anisotropy 
line 7 = (blue). Overlaid is a Gaussian distribution func- 
tion with the same mean and variance (red) . (b) Distribution 
of In X for the same parameters (blue) along with the corre- 
sponding distribution function of the logarithm of a Gaussian 
random variable. Here a = 0.3, L = 500, and A = 0.5. The 
distribution is non-Gaussian with a power-law tail towards 
large fidelity susceptibilities. 
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FIG. 14: (color online) (a) Distribution for 7 = 0.03 (near the 
value of 7 coinciding with the peak in [xli^p) (blue). Overlaid 
is a Gaussian distribution function with the same mean and 
variance (red), (b) Distribution of lux for the same parame- 
ters (blue) along with the corresponding distribution function 
of the logarithm of a Gaussian random variable. Here a = 0.3, 
L — 500, and A = 0.5. Note that there does not appear to be 
a power-law tail towards large fidelity susceptibilities. 



IV. CONCLUSION 

In this work we have studied the effect of random trans- 
verse fields and couplings on the phase diagram of the 
quantum XY chain. By examining the finite-size scaling 
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FIG. 15: (color online) (a) Distribution P(lnx)) at 7 = for 
L G {200,300,400,500}. (b) Rescaled distribution, ln(x/i) 
for the same parameters. Here a = 0.3 and A = 0.5, though 
this collapse applies for all other disorder strengths consid- 
ered. 
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FIG. 16: (color online) (a) Distribution P(lnx)) at 7 = 0.06, 
the position of the maximum of [xlt^p for = 0.3. System 
sizes L e {200, 300, 400, 500} are plotted, (b) Rescaled dis- 
tribution, ln(L^°x) for the same parameters, with a = 1.45. 
Here a — 0.3 and A = 0.5. 



tailed results suggest that the fidehty susceptibihty may 
be a useful tool for the study of other disordered systems. 
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V. APPENDIX 

Here we derive an expression for the fidelity suscepti- 
bility X in terms of the unitary matrix T: 



F{Z, Z) = \ \ dct 



T + T, 



exp < Tr In 



1 +Ttf 



1/2- 



exp Trln(i±llf±i^) 



1/2' 



exp { Tr In 1 



exp ■|^'^'^2 ^ ^ 



2 



1/2- 



exp { Tr- 



2 8 



(7) 



where 6T = dn-Tdx. This leads to 



of the typical fidelity susceptibility and the fidelity sus- 
ceptibility distribution for a range of disorder strengths 
and system sizes, we find agreement with earlier analytic 
results pertaining to the limit of weak randomness. The 
introduction of disorder clearly removes the anisotropy 
quantum critical line, replacing it with an extended Grif- 
fiths phase. There, the typical fidelity susceptibihty's 
finite-size scaling dimension depends strongly on the av- 
erage value of the anisotropy parameter, and appears to 
become extensive at the line of vanishing anisotropy. At 
the Ising critical line, the stretched exponential scaling of 
the fidelity susceptibility distribution is consistent with 
what is expected of an infinite randomness fixed point, 
while a Griffiths phase is observed to form in the vicin- 
ity. Remarkably, the scaling of the fidelity susceptibility 
distribution at the Ising critical line is universal, in that 
all disorder strengths give the same scaling behavior. In 
the Griffiths phase the fidelity susceptibility is not self- 
averaging. However, self-averaging behavior returns suf- 
ficiently far from the disorder-free critical lines. These de- 



F(Z,Z)=exp|Tri 



4 ^ 8 ^ ^ 



= exp 
= exp 
— exp 
= exp 



\d^T\\ldx'' 



1 



d^T\\%dx^ 



d.^T\\ldx^ 



\d^T\\ldx^ 



1 

16 
1 

16 
1 

16 



Tr [T^fd-^TT^d^T] dx^ 
Tr [-d:,T''TT''d^T] dx^ 
Tr [-d^T^'d^T] dx^ 



(8) 



where we have used the anti-symmetry of T^dxT which 
implies Tr [T-^a^T] = -||a^T|||,, with || • \\f the Frobe- 
nius norm. From ^ it follows that 



X=z\\9.T\\f- 



(9) 
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